Rejection-free Monte Carlo sampling for general potentials.
نویسندگان
چکیده
A Monte Carlo method to sample the classical configurational canonical ensemble is introduced. In contrast to the Metropolis algorithm, where trial moves can be rejected, in this approach collisions take place. The implementation is event-driven; i.e., at scheduled times the collisions occur. A unique feature of the new method is that smooth potentials (instead of only step-wise changing ones) can be used. In addition to an event-driven approach, where all particles move simultaneously, we introduce a straight event-chain implementation. As proof of principle, a system of Lennard-Jones particles is simulated.
منابع مشابه
Monte Carlo Techniques for Bayesian Statistical Inference – A comparative review
In this article, we summariseMonte Carlo simulationmethods commonly used in Bayesian statistical computing. We give descriptions for each algorithm and provide R codes for their implementation via a simple 2-dimensional example. We compare the relative merits of these methods qualitatively by considering their general user-friendliness, and numerically in terms of mean squared error and computa...
متن کاملMonte Carlo Integration With Acceptance-Rejection
This article considers Monte Carlo integration under rejection sampling or Metropolis-Hastings sampling. Each algorithm involves accepting or rejecting observations from proposal distributions other than a target distribution. While taking a likelihood approach, we basically treat the sampling scheme as a random design, and define a stratified estimator of the baseline measure. We establish tha...
متن کاملRejection-free kinetic Monte Carlo simulation of multivalent biomolecular interactions
The system-level dynamics of multivalent biomolecular interactions can be simulated using a rule-based kinetic Monte Carlo method in which a rejection sampling strategy is used to generate reaction events. This method becomes inefficient when simulating aggregation processes with large biomolecular complexes. Here, we present a rejectionfree method for determining the kinetics of multivalent bi...
متن کاملOn sequential Monte Carlo, partial rejection control and approximate Bayesian computation
We present a sequential Monte Carlo sampler variant of the partial rejection control algorithm introduced by Liu (2001), termed SMC sampler PRC, and show that this variant can be considered under the same framework of the sequential Monte Carlo sampler of Del Moral et al. (2006). We make connections with existing algorithms and theoretical results, and extend some theoretical results to the SMC...
متن کاملAn iterative Monte Carlo method for nonconjugate Bayesian analysis
The Gibbs sampler has been proposed as a general method for Bayesian calculation in Gelfand and Smith (1990). However, the predominance of experience to date resides in applications assuming conjugacy where implementation is reasonably straightforward. This paper describes a tailored approximate rejection method approach for implementation of the Gibbs sampler when nonconjugate structure is pre...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Physical review. E, Statistical, nonlinear, and soft matter physics
دوره 85 2 Pt 2 شماره
صفحات -
تاریخ انتشار 2012